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Energy  density  and  power  density  are  two  of  the  most  significant  performance  indices  of  a  fuel  cell  sys¬ 
tem.  Both  the  indices  are  closely  related  to  the  operating  conditions.  Energy  density,  which  can  be  derived 
from  fuel  cell  efficiency,  is  especially  important  to  small  and  portable  applications.  Generally  speaking, 
power  density  can  be  easily  obtained  by  acquiring  the  voltage  and  current  density  of  an  operating  fuel 
cell.  However,  for  a  direct  methanol  fuel  cell  (DMFC),  it  is  much  more  difficult  to  evaluate  its  efficiency 
due  to  fuel  crossover  and  the  complex  architecture  of  fuel  circulation.  The  present  paper  proposes  a 
semi-empirical  model  for  the  efficiency  evaluation  of  a  DMFC  under  various  operating  conditions.  The 
power  density  and  the  efficiency  of  a  DMFC  are  depicted  by  explicit  functions  of  operating  temperature, 
fuel  concentration  and  current  density.  It  provides  a  good  prediction  and  a  clear  insight  into  the  relation¬ 
ship  between  the  aforementioned  performance  indices  and  operating  variables.  Therefore,  information 
including  power  density,  efficiency,  as  well  as  remaining  run-time  about  the  status  of  an  operating  DMFC 
can  be  in  situ  evaluated  and  predicted.  The  resulting  model  can  also  serve  as  an  important  basis  for 
developing  real-time  control  strategies  of  a  DMFC  system. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Energy  density  and  power  density  are  two  of  the  most  sig¬ 
nificant  specifications  of  a  fuel  cell  system.  These  performance 
indices  usually  correlate  with  commercial  benefits  of  the  products 
[1,2].  Generally  speaking,  power  density  can  be  easily  obtained  by 
acquiring  the  voltage  and  current  density  of  an  operating  fuel  cell. 
However,  for  a  DMFC,  it  can  be  much  more  complex  to  evaluate  its 
energy  density  or  alternatively  its  efficiency.  The  real  efficiency  of  a 
fuel  cell  is  always  less  than  the  reversible  thermodynamic  efficiency 
due  to  voltage  losses  and  fuel  utilization  losses.  When  hydrogen  fuel 
is  considered,  the  efficiency  of  a  fuel  cell  can  be  evaluated  simply 
based  on  its  operating  voltage,  generated  current  and  the  flow  rate 
of  fuel  [3].  An  overabundance  of  fuel  supplying  to  a  fuel  cell  with¬ 
out  fuel  recycling  can  lead  to  a  waste.  Nonetheless,  for  a  DMFC,  the 
architecture  of  fuel  circulation  is  quite  complicated,  which  makes 
difficulties  in  evaluating  fuel  utilization  losses.  Fuel  crossover  is  one 
of  the  issues  reducing  the  voltage  efficiency  and  the  fuel  utilization 
efficiency  of  a  DMFC  [4].  To  suppress  the  fuel  crossover,  the  fuel  con¬ 
centration  of  a  DMFC  should  be  appropriately  regulated  between 
1.5  M  and  2.0  M  to  achieve  optimum  fuel  cell  efficiency  [5].  Taking 
the  energy  density  into  account,  fuel  of  higher  concentration  is  usu¬ 
ally  stored  in  a  fuel  cell  system.  The  fuel  is  then  diluted  in  a  mixing 
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reservoir  before  being  delivered  into  the  fuel  cell  stack.  The  resid¬ 
ual  (unreacted)  liquid  fuel  flowing  through  the  anode  is  circulated 
back  to  the  mixing  reservoir  and  then  to  the  stack  again.  Accord¬ 
ingly,  the  fuel  utilization  efficiency  cannot  be  estimated  merely  by 
the  generated  current  and  the  fuel  flow  rate.  It  is  therefore  essential 
to  determine  the  amount  of  fuel  crossover  losses  in  a  fuel  cell  stack. 

The  efficiency  of  a  fuel  cell  is  closely  related  to  its  operating 
conditions  such  as  fuel  concentration,  operating  temperature,  cur¬ 
rent  density  (or  operating  voltage),  rates  of  fuel  delivered,  etc.  Jiang 
et  al.  proposed  a  method  of  fuel  circulation  with  a  fixed  amount  of 
fuel  to  investigate  the  efficiency  of  a  DMFC  under  various  operating 
temperatures  [6].  Moreover,  the  influences  of  the  initial  fuel  con¬ 
centration  and  the  operating  temperature  on  the  operating  voltage 
of  a  passive  (without  fuel  and  air  circulation)  DMFC  were  also  stud¬ 
ied  [7].  Liu  et  al.  investigated  the  effect  of  the  membrane  thickness 
on  the  efficiency  of  a  passive  DMFC  under  various  operating  con¬ 
ditions  [8].  Seo  and  Lee  experimentally  determined  the  methanol 
crossover  fluxes  and  the  efficiency  of  a  DMFC  where  the  cell  tem¬ 
perature,  methanol  concentration,  aqueous  fuel  flow  rate,  air  flow 
rate  and  cathode  backpressure  were  considered  [9].  The  aforemen¬ 
tioned  articles  were  mainly  based  on  experimental  approaches. 
Although  these  experimental  approaches  were  straightforward,  it 
was  time-consuming  and  rather  difficult  to  depict  a  whole  picture 
of  the  efficiency  under  various  operating  conditions. 

On  the  other  hand,  theoretical  models  are  much  more  exact  and 
can  describe  details  and  provide  thorough  information  concerning 
the  system  [10-12].  Most  of  them  are  capable  of  forming  a  firm 
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Nomenclature 

List  of  symbols 

c 

bulk  methanol  concentration  in  the  channel 
(mol  cm-3) 

Da 

effective  diffusivity  of  methanol  in  anode  catalyst 
layer  (cm2  s-1 ) 

Db 

effective  diffusivity  of  methanol  in  anode  backing 
layer  (cm2  s-1 ) 

Dm 

effective  diffusivity  of  methanol  in  membrane 
(cm2  s-1) 

Er 

thermodynamically  reversible  voltage  (V) 

F 

Faraday’s  constant  (C  mol-1 ) 

Agrxn 

the  change  of  the  Gibbs  free  energy  of  the  reaction 
(J  mol-1) 

Afarxn 

the  change  of  the  enthalpy  of  the  reaction  (J  mol-1 ) 

j 

current  density  (A  cm-2) 

I< i 

partition  coefficient  at  the  interface  between  anode 
backing  layer  and  anode  catalyst  layer 

K 2 

partition  coefficient  at  the  interface  between  anode 
catalyst  layer  and  membrane 

Neff 

molar  flux  of  the  effective  fuel  consumption 
(mol  cm-2  s-1) 

Nmcr 

methanol  crossover  flux  (mol  cm-2  s-1 ) 

nd 

electro-osmotic  drag  coefficient  of  water 

nf 

remaining  amounts  of  the  fuel  (mol) 

Pdmfc 

power  density  of  a  direct  methanol  fuel  cell 

(Won-2) 

Re 

area-specific  resistance  (£2  cm2) 

T 

temperature  (I<) 

tr 

remaining  run-time  (h) 

tr 

specific  remaining  run-time  (hem2  mol-1 ) 

V 

operating  fuel  cell  voltage  (V) 

8A 

thickness  of  anode  catalyst  layer  (cm) 

8b 

thickness  of  anode  backing  layer  (cm) 

8m 

thickness  of  membrane  (cm) 

^dmfc 

efficiency  of  a  direct  methanol  fuel  cell 

£fuel 

fuel  utilization  efficiency 

^thermo 

reversible  thermodynamic  efficiency 

^■voltage 

voltage  efficiency 

X 

water  content  in  Nation 

a 

ionic  conductivity  of  membrane  (£2  cm2) 

basis  to  investigate  the  efficiency  of  a  DMFC.  Nonetheless,  they  are 
inevitably  complex  and  require  obscure  numerical  schemes  for  the 
solutions.  Besides,  parameters  of  the  model  should  be  measured  or 
identified  beforehand.  Such  a  task  is  usually  not  easy;  meanwhile, 
further  experimental  validation  is  essential.  Due  to  the  complexi¬ 
ties  of  the  solution  process,  theoretical  models  may  be  less  practical 
for  real-time  applications.  Accordingly,  semi-empirical  approaches 
are  anticipated  to  correlate  the  experimental  studies  with  the  the¬ 
oretical  models  that  can  serve  as  the  basis  for  developing  control 
strategies  of  a  fuel  cell  system. 

Since  the  voltage  losses  and  the  fuel  utilization  losses  are  the 
most  significant  factors  determining  the  efficiency  of  a  DMFC, 
we  briefly  review  the  polarization  and  the  fuel  crossover  semi- 
empirical  models.  The  semi-empirical  models  for  the  polarization 
behavior  of  a  DMFC  had  been  extensively  investigated.  Srini- 
vasan  et  al.  launched  the  issue  by  employing  activation  and  ohmic 
losses  to  depict  the  polarization  curves  of  a  polymer  electrolyte 
membrane  fuel  cell  (PEMFC)  [13].  Kim  et  al.  then  introduced  an 
additional  exponential  term  to  govern  the  mass  transport  phe¬ 
nomenon  and  had  a  better  fit  to  experimental  results  [14].  The 
mass  transport  term  was  modified  by  Squadrito  et  al.  to  provide 


a  limit  to  the  available  current  density  [15].  Moreover,  Kim’s  model 
was  also  modified  and  validated  on  a  30-cell  stack  by  Chu  and 
Jiang  in  which  the  physical  meaning  of  the  parameters  was  fur¬ 
ther  discussed  [16].  Because  the  mass  transport  limitations  of  a 
DMFC  predominantly  occur  at  the  anode  which  differ  from  the 
case  of  a  hydrogen  fuel  cell,  Scott  et  al.  studied  the  effects  of 
methanol  concentration  on  the  anode  polarization,  together  with 
an  empirical  open  circuit  voltage  model  and  a  cathode  overpoten¬ 
tial  model,  to  predict  the  overall  cell  voltage  [17].  Scott’s  group 
then  introduced  Kim’s  and  Squadrito’s  models  to  a  DMFC  [18]. 
It  revealed  that  Squadrito’s  model  fitted  the  results  better  when 
the  methanol  concentration  was  lower  than  0.75  M.  Furthermore, 
a  simple  equation  based  on  Squadrito’s  model  was  proposed  and 
discussed  by  Scott’s  group  [19].  To  elucidate  the  effect  of  methanol 
crossover  on  a  DMFC,  O’Hayre  et  al.  introduced  an  additional  term 
to  Scott’s  model  to  describe  voltage  losses  due  to  current  leakage 
and  fuel  crossover  [3].  Tu  et  al.  proposed  a  more  complicated  semi- 
empirical  model  for  distinguishing  the  individual  voltage  losses  due 
to  methanol  crossover  and  the  cathode  and  anode  overpotentials 
[20].  All  the  aforementioned  semi-empirical  models  incorporated 
certain  parameters  whose  values  could  be  influenced  by  operat¬ 
ing  variables  like  temperatures  and  fuel  concentrations.  Flowever, 
those  parameters  had  not  been  expressed  as  functions  of  operating 
variables  in  the  fitting  procedures.  Therefore,  they  are  not  suitable 
for  predicting  polarization  behavior  within  a  multi-dimensional 
domain  of  various  operating  variables. 

The  other  essentiality  of  the  efficiency  model  is  the  evaluation 
of  fuel  crossover.  Semi-empirical  models  for  evaluating  methanol 
crossover  fluxes  had  been  extensively  discussed  [21-24].  In  our 
previous  work,  a  methanol  crossover  model  was  established  ana¬ 
lytically  using  an  algebraic  function  in  which  fuel  concentration, 
current  density  and  operating  temperature  of  a  DMFC  were  con¬ 
sidered  [25].  It  provides  a  good  prediction  and  a  clear  picture  of 
the  methanol  crossover  behavior  in  a  three-dimensional  domain  of 
operating  variables.  Consequently,  the  model  is  employed  in  this 
study  for  efficiency  evaluation. 

The  major  objective  of  this  work  is  aimed  at  the  efficiency  evalu¬ 
ation  of  a  DMFC  in  a  multi-dimensional  domain  of  various  operating 
variables.  The  proposed  semi-empirical  model  takes  performance 
indices  (i.e.  power  density  and  efficiency)  and  operating  variables 
(i.e.  fuel  concentration,  current  density  and  temperature  of  the  fuel 
cell)  into  a  whole  picture  simultaneously.  In  the  model,  both  the 
power  density  and  the  efficiency  are  stated  in  analytical  functions 
which  explicitly  incorporate  fuel  concentration,  current  density 
and  temperature  of  the  fuel  cell.  Once  the  corresponding  param¬ 
eters  of  the  model  are  determined,  the  quantitative  prediction  of 
power  density  and  efficiency  can  be  easily  obtained.  As  a  result,  it 
can  serve  as  an  important  basis  for  developing  optimum  controlling 
strategies  of  a  DMFC  [26].  The  model  can  also  predict  the  remain¬ 
ing  run-time  of  a  DMFC  system  under  various  operating  conditions. 
Such  information  is  also  useful  for  practical  situations. 


2.  Model  description 

The  behavior  of  a  fuel  cell  is  closely  related  to  its  operating 
conditions,  i.e.  operating  voltage  (or  current),  operating  tempera¬ 
ture,  fuel  concentration,  flow  rates  of  fuel  and  air,  and  the  pressure 
built  in  the  anode  and  cathode  chambers.  In  this  paper,  fuel  con¬ 
centration  c,  current  density  j,  and  operating  temperature  T  are 
considered  as  the  significant  factors  of  the  semi-empirical  model 
for  evaluating  the  DMFC  efficiency  £dmfc  under  steady-state  condi¬ 
tions.  In  the  following  sections,  preliminaries  of  the  definitions  and 
assumptions  of  £dmfc  will  be  addressed  first.  Then  a  semi-empirical 
model  of  polarization  behavior  is  proposed.  By  incorporating  an 
algebraic  semi-empirical  model  for  fuel  crossover  evaluation,  the 
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relationship  between  power  density  Pdmfc  and  efficiency  £dmfc  and  Vofa  fuel  cell.  The  present  work  is  based  on  the  model  proposed 

under  various  operating  conditions  is  established.  by  Scott’s  group  [18,19]: 


2.1.  Preliminaries 

The  £dmfc  mainly  consists  of  reversible  thermodynamic  effi¬ 
ciency  £  thermo*  voltage  efficiency  of  the  fuel  cell  £voitage*  and  fuel 
utilization  efficiency  £fuei  [3].  It  can  be  written  as: 


^dmfc  —  ^thermo  x  ^-voltage  x  ^fuel 


(1) 


First,  the  reversible  thermodynamic  efficiency  is  defined  as: 


0  _  Agrxn 

^thermo  — 

zAUrxn 


(2) 


where  Agrxn  and  Ahrxn  are  the  changes  of  the  Gibbs  free  energy 
and  the  enthalpy,  respectively,  of  the  reaction.  When  the  assump¬ 
tion  of  an  isobaric  condition  during  the  reaction  is  made,  both  the 
Gibbs  free  energy  and  the  enthalpy  of  reaction  are  temperature 
dependent.  Considering  the  oxidation  reaction  of  a  DMFC  at  anode: 

CH3OH  +  H20  -*  C02  +  6H+  +  6e"  (3) 


and  the  reduction  reaction  at  cathode: 


§02  +  6H+  +  6e-  -*  3H20  (4) 

the  £thermo  of  a  DMFC  is  around  96.7-96.6%  when  T  is  within  the 
range  of  298  K  to  340  K.  For  the  reason  of  simplicity,  we  define 
£ thermo  =  96.7  0/0  and  neglect  the  temperature  effect  in  the  following 
illustration  [3,27]. 

As  for  the  voltage  efficiency  £voitage  of  a  fuel  cell,  it  is  defined  as: 


_  V 

^■voltage  — 


(5) 


where  V  and  Er  denote  the  operating  voltage  and  the  thermody¬ 
namically  reversible  voltage,  respectively,  of  a  fuel  cell.  The  effects 
of  T  and  c  on  Er  can  be  evaluated  using  Nernst  equation.  It  indicates 
that  the  Er  of  a  DMFC  varies  between  about  1.20  V  and  1.21  V  at  a 
standard  pressure  (1  atm)  and  a  temperature  between  298  K  and 
340  K,  where  the  c  is  ranging  from  0.75  M  to  1.5  M.  Both  a  higher 
c  and  a  lower  T  will  lead  to  a  higher  Er.  In  this  work,  Er  =  1.21  V  is 
employed  under  a  feasible  operating  condition.  One  can  modify  Er 
as  a  function  of  T  and  c  by  introducing  Nernst  equation  [3].  More¬ 
over,  the  operating  voltage  V  will  be  interpreted  as  a  function  of  T, 
c  andj  in  Section  2.2. 

Finally,  the  £fuei  accounts  the  amounts  of  fuel  that  effectively 
generate  electricity  in  the  electrochemical  reaction  and  the  losses 
wasted  due  to  fuel  crossover.  As  a  result,  £fuel  is  expressed  as: 


V  =  Er  +  b  log  jo  -  b  log  j  -  Rej  +  Q  ln(l  -  CJ)  (8) 


where  b ,  j0  and  Re  are  the  Tafel  slope,  the  exchange  current 
density  and  the  area-specific  resistance,  in  order.  Meanwhile, 
Ci  and  C2  are  empirically  determined  parameters  regarding  the 
mass  transfer  overpotential.  The  above  model  had  been  verified 
and  was  found  to  fit  the  behavior  of  PEMFCs  and  DMFCs  very 
well.  However,  all  the  aforementioned  parameters  are  T  and  c 
dependent  [19].  The  form  of  Eq.  (8)  itself  cannot  be  used  to  pre¬ 
dict  the  polarization  behavior  under  various  operating  conditions, 
whereas  it  is  essential  for  developing  the  controlling  strategies  of 
a  fuel  cell  system.  In  this  study,  an  additional  term  r]cr  regard¬ 
ing  the  overpotential  due  to  fuel  crossover  is  added  to  Eq.  (8). 
The  relevant  parameters  log/o,  b,  Ci,  C2  and  rjcr  are  empirically 
approximated  by  cubic  polynomials  of  T  and  c.  Moreover,  Re  is 
expressed  as  a  temperature  function.  It  is  found  that  identifying  Re 
separately  is  helpful  to  enhance  numerical  stability  and  accuracy 
for  the  fitting  procedures  when  various  operating  conditions  are 
considered. 

Since  the  V  is  generally  measured  by  clamping  the  current  col¬ 
lectors  on  both  sides  of  the  fuel  cell  stack  (or  single-cell  module), 
the  Re  can  comprise  the  portions  regarding  ionic  charge  transport 
(mainly  within  the  electrolyte  membrane)  and  electronic  charge 
transport  (mainly  induced  between  the  interfaces  of  each  compo¬ 
nents  of  the  stack).  Nonetheless,  the  ionic  charge  transport  tends  to 
be  more  difficult  than  electronic  charge  transport  [3].  In  our  exper¬ 
iments  where  a  single-cell  DMFC  module  (35  mm  x  35  mm  MEA) 
is  adopted,  the  electronic  resistance  (without  MEA  and  gaskets  in 
the  testing  module)  was  less  than  7%  of  the  resistance  measured 
on  a  regular  operating  DMFC.  The  temperature  range  was  between 
303  K  and  333  K.  In  addition,  the  tendency  of  the  electronic  resis¬ 
tance  variation  against  the  temperature  is  similar  to  that  of  the 
ionic  resistance.  Accordingly,  the  Re  in  Eq.  (8)  is  approximated 
by  introducing  the  form  of  ionic  conductivity  a  of  the  electrolyte 
membrane  and  pursuing  a  best  fit  that  can  moderately  incorporate 
the  effect  of  electronic  charge  transport.  The  Re  can  be  expressed 
as: 

Re  =  ^  (9) 

<7 

where  Sm  is  the  thickness  of  the  electrolyte  membrane.  The 
ionic  specific  conductivity  of  Nafion-type  membrane  (duPont  de 
Nemours,  Wilmington,  DE)  can  be  expressed  as  a  function  of  T[28]: 


cr(T)  =  (0.5139A  -  0.326)  exp 


,26s  (as 


(10) 


£fuel 


Neff 

Neff  +  NmCr 


(6) 


where  Neff  and  Nmcr  represent  the  molar  fluxes  of  the  effective  fuel 
consumption  and  the  fuel  crossover  losses,  respectively.  And  Neff 
can  be  addressed  as  a  function  of  current  density  j.: 


Neff  =5=  (7) 

where  z  =  6  for  a  DMFC,  and  F  is  the  Faraday’s  constant.  Moreover, 
the  fuel  crossover  fluxes  had  been  investigated  in  our  previous  work 
[25].  It  will  be  briefly  reviewed  in  Section  2.3. 


where  X  is  the  water  content  of  a  membrane.  When  the  mem¬ 
brane  is  well-hydrated,  the  water  content  is  constant  and  the  ionic 
conductivity  can  be  defined  by  the  form  of  Eq.  (10): 

Because  8m  is  known  and  Re  can  be  easily  acquired  while  performing 
j-V measurements  under  various  T r\  and  r2  can  then  be  determined 
using  a  nonlinear  least-squares  algorithm  described  in  Appendix  A. 
Consequently,  Re{T)  can  be  interpreted  as: 


Re(T)  =  ^  exp  (y 


1268  \ 
303  ) 


(12) 


2.2.  Polarization  model 

Due  to  a  variety  of  overpotentials,  the  actual  operating  voltage  V 
of  a  fuel  cell  will  always  be  lower  than  its  Er.  There  are  various  semi- 
empirical  models  in  literature  to  depict  the  relationship  between  j 


The  operating  voltage  of  a  DMFC  can  be  influenced  by  the  cathode 
mixed  potential  due  to  methanol  crossover.  It  had  been  verified 
that  the  influence  of  methanol  on  the  cathode  mixed  potential  is 
independent  of  j  [17].  Therefore,  the  r\CT  as  well  as  b,  j0,  C\  and  C2 
are  interpreted  as  functions  of  T and  cThen  Eq.  (8)  can  be  expressed 
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as: 

V(T,  cj )  =  Er  +  b(T ,  c)  •  [a(T,  c)  -  logj]  -  Re(T)  j 

+C\(T,  c)  •  In  [1  -  C2(T,  c)  •  j]  -  Tjcr(T,  c)  (13) 

where  a(T,  c)  is  substituted  for  log[/o(T,  c)]  to  enhance  numeri¬ 
cal  robustness  of  the  fitting  algorithm.  The  functions  a(T,c),  b(T,c), 
Ci(T,c),  C2(T,c)  and  rjcr(T,  c)  are  empirically  approximated  by  cubic 
polynomials: 


Xj(T,  c)  =  +  a2iT  +  a3jC  +  a^T2  +  a5i7c  +  a6ic2  +  a7iT 3 

+a8iT2c  +  otQjTc2  +  a10iC3  (14) 

wherexj(i  =  l,2, . .  .,5)  can  denote  a,b,Ci,C2  and  r\cr-  Each  of  the  five 
functions  consists  of  10  parameters  to  be  fitted.  For  a  pre-assigned 
combination  of  (Tk,  ck),  Eq.  (13)  can  be  expressed  as: 

V(Tk,  ck,j)  =  Er  +  bk ■  [ak  -  log  j]  -  Re(Tk )  j  +  Cljk  •  In  [l  -  C2,k  •  j] 

—  bcr,k  (15) 


According  to  the  model,  the  functions  DA(T),  DB(T),  Dm(T)  and 
nd(T)  are  defined  as: 


Da(T)  =  fa  exp  fa  (3^3  -  ^) 
DB(T)  =  fa  x  10 -l-4163-/34/T 

Dm(T)  =  fa  exp  /36  (^3  - 


and 

nd(T)  =  2.9  exp 


Ms-f) 


(19) 

(20) 
(21) 

(22) 


In  Eqs.  (19)-(22),  there  are  eight  unknown  parameters  fy  (i  =  l, 
2,  . . .,  8)  to  be  fitted  by  adopting  the  Levenberg-Marquardt  algo¬ 
rithm  described  in  Appendix  A.  Once  the  eight  fitted  parameters 
are  determined,  the  quantitative  prediction  of  the  Nmcr(T,  c,  j)  will 
be  straightforward  with  no  need  of  any  numerical  strategy. 


2.4.  Power-efficiency  relationship 


The  j-V data  corresponding  to  each  (Tk,  ck)  acquired  experimentally 
was  substituted  into  Eq.  (15)  to  obtain  parameters  ak ,  bk,  clk,  c2k 
and  r)cr  k  using  Levenberg-Marquardt  algorithm,  which  was  briefly 
illustrated  in  Appendix  A.  Based  on  a  number  of  polarization  mea¬ 
surements  (for  k=  1,2,...,  n),  the  fitted  results  of  Eq.  (15)  together 
with  the  corresponding  (Tk,  ck)  are  then  substituted  into  Eq.  (14) 
to  determine  the  unknown  parameters  a^-,  a2i,  . . .,  a10l  using  the 
ordinary  least-squares  method  in  Appendix  B.  Consequently,  the 
operating  voltage  V  can  then  be  expressed  by  an  explicit  function 
of  T,  c  and  j,  as  was  shown  in  Eqs.  (12)-(14). 

2.3.  Methanol  crossover  model 

Fuel  crossover  phenomenon  is  considered  one  of  the  most  sig¬ 
nificant  factors  contributing  to  £fuel  of  a  DMFC.  In  our  previous 
work  [25],  the  methanol  crossover  flux  was  stated  analytically 
in  an  algebraic  function,  which  explicitly  incorporated  c,  j  and  T 
of  the  fuel  cell.  The  model  took  various  physical  parameters  into 
account  including  the  thicknesses  of  membrane,  catalyst  layer  and 
gas  diffusion  layer  (<$m,  8A,  and  <5B),  as  well  as  electro-osmotic  drag 
coefficient  of  methanol  nd  and  effective  diffusivities  ( Dm,DA  and  DB) 
at  the  corresponding  positions  in  a  membrane  electrode  assembly 
(MEA).  The  transfer  behavior  of  methanol  was  governed  by  diffu¬ 
sion  and  electro-osmosis.  The  aforementioned  coefficients  Dm,DA, 
Db  and  nd  were  also  expressed  by  functions  of  T  and  c.  By  adopt¬ 
ing  this  methanol  crossover  model,  the  molar  flux  of  the  methanol 
crossover  losses  Nmcr  in  Eq.  (6)  can  be  expressed  as: 

exp  \j8mK{T)  1 

Nmcr(T,  c,  j)  =  -j  ■  c2  m(c,j,  T)  ■  Dm(T )  •  k(T)- - 1  r.c  *  (16) 

1  -  exp  \j8mK{T)\ 


where 


Based  on  the  aforementioned  illustration,  the  £dmfc  can  then  be 
expressed  by  an  explicit  function  of  operating  variables  T,  c  and  j: 


£dmfc(T,  C,  j)  =  96.7%  x  V(J\C:j)  x 


J/6F 


(23) 


1.21  Q-/6  F)-Nmcr(T,c,j) 

In  addition,  the  power  density  of  a  DMFC  can  be  expressed  as: 
Pdmfc(T,cJ)=j  V{T,c,j)  (24) 


By  incorporating  Eqs.  (23)  and  (24),  a  clear  picture  that  simultane¬ 
ously  depicts  the  relationships  between  the  performance  indices 
(£dmfc  and  Pdmfc)  and  the  operating  variables  (T,  c  and  j)  is  revealed. 

Furthermore,  the  remaining  run-time  of  a  DMFC  can  be  pre¬ 
dicted  from  £dmfc  and  pdmfc  by: 


,  fT  n  =  nf  - Ah  rxn  '  £dmfc(^  C,  j) 

K  Pdmfc(F,  C,j)-A  -  3600 


(25) 


where  and  A  are  the  remaining  amount  of  fuel  in  the  fuel  tank  and 
the  effective  area  of  the  MEA,  respectively.  It  can  be  alternatively 
defined  based  on  unit  mole  of  the  fuel  and  the  specific  area  of  the 
MEA  as: 


_  A/lrxn  •  £dmfc(l^  C,  j) 
Pdmfc(^  C?j)  '  3600 


(26) 


Such  information  varies  with  operating  conditions  and  is  generally 
useful  for  practical  situations. 


3.  Experimental 

The  experiment  of  the  present  work  consists  of  two  parts.  One 
is  the  j-V  measurements  for  fitting  and  validation  of  the  pro¬ 
posed  polarization  model.  The  other  is  the  methanol  crossover 


K2  {-1  +exp  [j8mK(T)]  }  {Db(T)  [\2FcfaDA(T)-j8A]  -2 ]fa8BDA(T)} 

12 F  {— D/\(T)Dg(T)  +  exp  [jSmK(T)}  {Da(T)Db(T)  +jK2Dm(T)K(T )  [5ADB(T)  +  MBDA(T)]  }  } 


(17) 


is  the  methanol  concentration  at  the  interface  between  the  anode 
catalyst  layer  and  the  membrane,  K\  and  I<2  are  partition  coeffi¬ 
cients,  while  effective  diffusivities  DA,  DB  and  Dm  are  all  functions 
of  temperature  T.  Moreover 


measurements  for  establishing  the  semi-empirical  methanol 
crossover  model  and  thereby  the  relationship  between  £dmfC  and 
Pdmfc  defined  by  Eqs.  (23)  and  (24). 


nd(T)-s 
F>m(T)  ■  F 


(18)  3.1.  Experimental  facilities 


is  defined  for  the  simplicity  of  the  equation  expression,  where  the 
concentration  translation  factor  s  =  18.4  is  employed  for  diluted 
methanol  aqueous  solution. 


The  experimental  facilities  employed  in  this  study  were  the 
same  as  those  described  in  [25].  A  standard  testing  module  of 
a  unit  commercial  MEA  (35  mm  x  35  mm,  duPont  de  Nemours, 
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Table  1 

A  brief  list  of  the  experimental  conditions  for  the  fitting  and  the  validation  of  the  proposed  polarization  model. 
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Items 


Descriptions 


1. MEA 

2.  Graphite  channel  type 

3.  Operating  temperature  for  model  fitting 

4.  Fuel  concentration  for  model  fitting 

5.  Conditions  for  model  validation 

6.  Electric  load 

7.  Fixed  gas  flow  rate  on  the  cathode  side 

8.  Fixed  fuel  flow  rate  on  the  anode  side 

9.  Fuel  delivery  type 


35  mm  x  35  mm  (DuPont  Nation®  117) 

Single  serpentine,  1  mm  in  width  and  depth 
303  K,  313  K,  323  K,  333  K 
0.75  M,  1.0  M,  1.25  M,  1.5  M 

(a)  313  K:  0.875  M,  1.125  M,  1.375  M 

(b)  1.0  M:  308  K,  318  K,  328  K 

Potentiostatic  mode  from  0.6  V  to  0.15  V  with  an  increment  of  -0.05  V  and  the 
duration  of  5  min  for  each  measurement 
Feeding  air  with  1 50  mL  min-1 
5  mLmin-1 
Single  pass 


Wilmington,  DE)  was  adopted  for  the  experiments.  The  MEA  was 
Nafion-117  based.  The  anode  catalyst  was  Pt-Ru/C  (Pt-Ru  content 
40wt%,  Pt/Ru  wt  ratio  was  1/1)  and  the  cathode  catalyst  was  Pt/C 
(Pt  content  40wt%).  The  Pt  loadings  at  anode  and  cathode  of  the 
MEA  were  both  2  mg  cm-2.  The  flow  channel  was  single  serpen¬ 
tine  and  1  mm  in  width  and  depth.  A  temperature  control  unit  was 
applied  to  provide  appropriate  temperature  conditions.  The  fuel 
delivery  was  driven  by  a  precision  liquid  pump  (Eyela  MP-1000, 
Tokyo,  Japan).  In  order  to  keep  a  constant  fuel  concentration,  the 
fuel  delivery  was  single  pass  without  circulating  back  to  the  fuel 
tank.  To  avoid  the  influence  of  uncertainties  and  variation  of  C02 
concentration  in  the  atmosphere,  the  air  fed  into  cathode  was  pro¬ 
vided  by  a  gas  cylinder.  The  exhausted  stream  at  the  cathode  outlet 
was  guided  into  a  gas  chromatography  (China  Chromatography 
GC3000,  Taipei,  Taiwan)  with  a  thermal  conductivity  detector  (TCD) 
for  C02  measurement.  A  source  meter  (Keithley  2440,  Cleveland, 
OH)  was  introduced  to  provide  designate  electric  loads.  In  addition, 
a  precision  internal  resistance  meter  (Hioki  3561,  Nagano,  Japan) 
was  adopted  to  measure  the  area-specific  resistance  of  the  fuel  cell 
by  clamping  the  current  collectors  on  both  sides  of  the  single-cell 
module. 

3.2.  j-V  data  acquisition  for  the  polarization  model 

Four  levels  of  T  and  c  were  assigned  for  j-V  measurements  as 
shown  in  Table  1.  There  were  16  combinations  of  (T^,  Cjf)  with  k=  1, 
2,  . . .,  16.  The  j-V  measurements  were  carried  out  by  the  poten¬ 
tiostatic  technique  with  voltages  ranged  from  0.6  V  to  0.15  V  and  a 
decrement  of  0.05  V.  The  duration  of  each  measurement  was  5  min. 
Thus,  10  data  points  for  each  j-V  measurement  were  obtained. 
The  flow  rates  of  air  and  fuel  were  150  mLmin-1  and  5  mLmin-1, 
respectively.  The  internal  resistance  corresponding  to  each  (Tk,  C/J 
was  also  acquired  for  the  fitting  of  the  polarization  curve.  The 
resulting  data  were  then  introduced  to  establish  the  polarization 
model  following  the  procedure  described  in  Section  2.2. 

To  validate  the  resulting  polarization  model,  additional  six 
j-V  measurements  were  carried  out  at  a  constant  T=313K  with 
c  =  0.875  M,  1 .1 25  M  and  1 .375  M,  and  at  a  constant  c  =  1 .0  M  with 
T=308  K,  318  K  and  328  K.  These  measured  data  were  compared 


with  the  estimated  results  using  the  resulting  polarization  model 
to  verify  its  feasibility. 

3.3.  C02  measurements  for  the  methanol  crossover  model 

Methanol  crossover  measurements  were  carried  out  using  the 
same  procedures  as  those  described  in  [25].  Although  most  of  the 
methanol  that  transmits  to  cathode  can  be  oxidized  to  C02,  there 
are  possibilities  that  some  methanol  may  remain  unreacted  and 
induce  certain  intermediates.  Many  researchers  assumed  that  the 
amount  of  the  unreacted  methanol  was  negligible  when  comparing 
with  the  overall  crossover  methanol  [24,29],  whereas  Eccarius  et  al. 
suggested  the  residual  methanol  cannot  be  omitted  in  some  cases 
[22].  Therefore,  the  effect  of  the  residual  methanol  on  the  crossover 
evaluation  is  still  an  open  issue.  Moreover,  the  induced  intermedi¬ 
ates  are  generally  neglected,  since  the  amounts  are  very  small  and 
the  reaction  rates  are  fast  [30].  In  our  previous  work  [25],  the  influ¬ 
ence  of  CO  on  the  methanol  crossover  measurement  was  verified  to 
be  less  than  1%.  Accordingly,  it  is  assumed  in  this  paper  that  all  the 
methanol  that  transmits  to  cathode  is  instantaneously  oxidized  to 
C02  and  H02,  and  the  intermediates  such  as  CO,  CH2OH  and  CHOH 
are  neglected.  The  methanol  crossover  fluxes  can  be  quantitatively 
identified  by  measuring  the  C02  concentration  of  the  cathode  out¬ 
let  exhausted  stream.  A  brief  list  of  the  experimental  conditions  for 
establishing  the  methanol  crossover  model  was  shown  in  Table  2. 
Two  levels  of  T  and  c  were  assigned,  thus  obtaining  four  sets  of  T 
and  c  combinations.  For  each  set  of  T  and  c,  7  levels  ofj  equally  dis¬ 
tributed  within  the  feasible  region  were  assigned.  The  electric  load 
was  applied  by  the  galvanostatic  mode.  Accordingly,  there  were  28 
sets  of  C02  molar  fluxes  to  be  acquired  for  an  operating  DMFC. 

Nonetheless,  not  only  the  methanol  but  also  the  C02  produced 
at  the  anode  can  permeate  across  the  membrane  to  cathode.  Such 
amount  of  C02  should  be  subtracted  from  the  C02  measurement 
results  in  an  operating  DMFC  to  identify  the  methanol  crossover 
fluxes  more  accurately.  The  half-cell  measurement  scheme  was 
then  introduced  to  measure  the  C02  crossover  fluxes  [22,25,31], 
where  the  air  flow  on  the  cathode  side  was  alternatively  replaced 
by  an  inert  gas  such  as  nitrogen,  and  the  specified  currents  were 
drawn  by  a  source  meter.  In  such  situations,  even  some  methanol 


Table  2 

A  brief  list  of  the  experimental  conditions  for  the  fitting  and  the  validation  of  the  crossover  model. 

Items 

Descriptions 

1.  Operating  temperature  for  model  fitting 

2.  Fuel  concentration  for  model  fitting 

3.  Conditions  for  model  validation 

4.  Electric  load 

5.  Fixed  gas  flow  rate  on  the  cathode  side 

6.  Fixed  fuel  flow  rate  on  the  anode  side 

7.  Fuel  delivery  type 

303  K,  333  K 

1.0M,  2.0  M 

1.5  M:  303  K,  313  K,  323  K,  333  K 

Assigning  7  levels  of  current  density  equally  distributed  within  the  feasible 
region  corresponding  to  each  combination  of  temperature  and  concentration 

150 mLmin-1  (feeding  N2  for  half-cell  measurement,  air  for  ordinary  DMFC 
measurement) 

5  mLmin-1 

Single  pass 
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Temperature  /  K 

Fig.  1.  The  area-specific  resistances:  (a)  the  measured  values  under  various  levels  of 
fuel  concentration  and  operating  temperature:  (O)  303  K,  (A)  313  K,  (□)  323  K  and 
(v)  333  K;  (b)  the  measured  values  (symbols)  and  the  fitted  results  (curve). 

transmitted  through  the  membrane  from  anode  to  cathode,  it  did 
not  perform  any  reaction  due  to  the  lack  of  oxygen.  Therefore, 
no  any  additional  C02  was  produced  at  cathode.  If  any  C02  was 
detected  at  the  cathode  outlet  under  such  operating  modes,  it 
was  recognized  as  that  transmitted  from  anode,  thereby  verifying 
the  C02  crossover  fluxes.  The  C02  crossover  fluxes  in  the  half-cell 
mode  were  also  measured  according  to  the  same  conditions  listed 
in  Table  2.  Accordingly,  the  methanol  crossover  fluxes  could  be 
obtained  for  fitting  the  eight  unknown  parameters  /32-  (z  =  1,  2, 

8). 

It  is  noted  that  the  C02  crossover  fluxes  in  half-cell  measure¬ 
ments  may  not  be  identical  to  that  in  an  ordinary  operating  DMFC, 
thus  influencing  the  compensation  of  the  methanol  crossover  mea¬ 
surement.  The  C02  concentrations  at  the  anodes  of  both  cases  are 
the  same,  whereas  the  C02  concentration  at  the  cathode  under 
the  half-cell  condition  may  be  lower  than  that  under  the  ordi¬ 
nary  operating  condition.  It  is  because  C02  can  be  generated  at  the 
cathode  under  the  ordinary  operating  condition  only.  Fortunately, 
such  a  deviation  between  the  cathode  C02  concentrations  of  both 
cases  is  negligible  when  the  methanol  crossover  is  slight.  When 
the  methanol  crossover  is  severer,  our  previous  work  has  verified 
that  the  C02  crossover  becomes  less  significant  when  comparing 
with  the  methanol  crossover  [25].  Accordingly,  the  difference  of 
C02  crossover  between  half-cell  and  ordinary  operating  conditions 
is  neglected  in  this  work,  as  was  assumed  in  the  literatures  [22,31  ]. 

Finally,  to  validate  the  proposed  efficiency  model,  four  sets  of 
additional  methanol  crossover  measurements  were  carried  out 
at  a  constant  c=1.5M  with  T=303K,  313 K,  323 K  and  333 K.  For 
each  measurement  condition,  the  corresponding  j-V  measure¬ 
ments  were  also  performed  to  determine  £dmfc  and  Pdmfc  for  the 
comparison  with  the  estimated  results  using  the  achieved  effi¬ 
ciency  model. 

4.  Results  and  discussions 

41.  Polarization  model 

By  following  the  procedure  described  in  Section  3.2,  a  total  of 
160  j-V  data  were  acquired.  For  each  combination  of  ( Tk ,  ck),  the 
internal  resistance  was  also  obtained  and  shown  in  Fig.  1(a).  It 
revealed  that  the  internal  resistance  of  an  operating  DMFC  was 
independent  of  c,  whereas  higher  T  led  to  lower  internal  resis¬ 
tances.  The  results  shown  in  Fig.  1(a)  were  inserted  into  Eq.  (12) 


Concentration  /  M 

Fig.  2.  The  fitted  values  of  the  parameters  in  Eq.  (15)  corresponding  to  various  levels 
of  methanol  concentration  c  and  operating  temperature:  (O)  303  K,  (A)  313  K,  (□) 
323 1<  and  (v)  333 1<. 


to  obtain  the  parameters  r\  =  1.46  x  10-2  cm-1  and  r2  =  656.41 
K  as  Sm  =  0.018  cm.  The  experimental  data  and  the  fitting  results 
of  Eq.  (12)  are  shown  in  Fig.  1(b).  The  relative  error  between  the 
experimental  data  and  the  estimated  values  was  around  1.3%.  Such 
a  temperature  dependent  function  Re(T )  was  then  employed  in  Eqs. 
(13)  and  (15). 

For  each  combination  of  ( Tk ,  ck ),  the  corresponding  10  j-V  data 
were  inserted  into  Eq.  (15)  to  obtain  the  parameters  ak,  bk ,  Clfc,  C2k 
and  ricr  k  for  k=  1,  2, . . .,  16.  The  fitting  results  of  these  parameters 
are  illustrated  in  Fig.  2.  All  these  values  together  with  each  individ¬ 
ual  combination  of  ( Tk ,  ck)  were  introduced  to  Eq.  (14)  to  obtain  the 
parameters  al2,  a2l-, . . .,  a10l-  (z  =  1 ,  2  , . . .,  5)  by  means  of  linear  least 
squares  method  described  in  Appendix  B.  The  form  of  these  func¬ 
tions  and  the  values  of  the  fitted  parameters  are  summarized  in 
Table  B.l.  Accordingly,  the  semi-empirical  polarization  model  of  a 
DMFC  under  various  operating  conditions  was  then  accomplished. 

The  resulting  model  was  validated  using  additional  six  sets  ofj- V 
measurements  obtained  according  to  the  T  and  c  listed  in  Table  1. 
Fig.  3  shows  the  polarization  curves  calculated  using  Eq.  (13)  under 
various  operating  temperatures  with  a  fuel  concentration  of  1.0  M. 
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Fig.  3.  Polarization  curves  based  on  experimental  results  (symbols)  and  the  pro¬ 
posed  model  (curves).  The  fuel  concentration  is  1.0  M.  The  operating  temperatures 
are:  (O)  303  K,  (•)  308  K,  (A)  313  K,  (a)  318  K,  (□)  323  K,  (■)  328  K  and  (v)  333  K. 

The  empty  symbols  represent  the  experimental  data  that  were 
employed  for  fitting  the  unknown  parameters,  while  the  filled  sym¬ 
bols  are  the  additional  data  for  model  validation.  It  can  be  observed 
that  the  prediction  of  the  proposed  model  meets  the  experimental 
data  quite  well.  The  relative  deviation  between  the  predicted  and 
the  measured  voltages  was  around  3.3%.  Similarly,  the  proposed 
model  was  also  validated  under  various  concentration  levels.  Fig.  4 
shows  the  polarization  curves  with  the  fuel  concentration  ranging 
from  0.75  M  to  1.5  M,  whereas  the  operating  temperature  was  kept 
at  313  K.  It  reveals  a  good  coincidence  between  the  predicted  and 
the  measured  results.  The  relative  deviation  between  the  predictive 
and  measured  data  was  around  2.1%.  Consequently,  the  feasibility 
of  the  proposed  polarization  model  can  be  verified. 

4.2.  Power-efficiency  relationship 

Using  the  semi-empirical  methanol  crossover  model  described 
in  Section  2.3  and  the  procedures  mentioned  in  Section  3.3,  an  illus¬ 
tration  of  Nmcr(T,  c,  j)  in  the  space  of  multi-operating  variables  was 
obtained  and  shown  in  Fig.  5.  The  V(T ,  c,  j)  obtained  in  Section  4.1 
and  Nmcr(T,  c,  j)  obtained  in  this  section  were  then  introduced  to 


Fig.  4.  Polarization  curves  based  on  experimental  results  (symbols)  and  the  pro¬ 
posed  model  (curves).  The  operating  temperature  is  313  K.  The  fuel  concentrations 
are:  (O)  0.75  M,  (•)  0.875  M,  (A)  1.0  M,  (a)  1.125  M,  (□)  1.25  M,  (■)  1.375  M,  and  (v) 
1.5  M. 


Eq.  (23)  to  calculate  £dmfc-  An  illustration  of  the  £dmfc(^  c,  j)  in 
the  space  of  multi-operating  variables  is  depicted  in  Fig.  6.  A  com¬ 
plex  and  highly  nonlinear  behavior  coupled  with  multi-variables 
is  observed.  It  therefore  reveals  the  essentiality  of  the  proposed 
model  that  provides  a  clear  picture  of  £dmfc  influenced  by  vari¬ 
ous  operating  variables.  A  special  situation  at  333  K  is  extracted 
from  Fig.  6  to  illustrate  the  influence  of  c  on  £dmfc  and  is  shown  in 
Fig.  7.  An  increasing  £dmfc  with  j  is  observed  in  the  lower  j  region. 
It  is  because  the  methanol  crossover  fluxes  are  attenuated  by  an 
increase  in  current  density,  thus  enhancing  the  fuel  utilization  effi¬ 
ciency.  On  the  other  hand,  a  decrease  of  £dmfC  with  j  is  observed 
in  the  higher  J  region.  It  is  because  the  operating  voltage  is  low¬ 
ered  when  approaching  the  limiting  current  density,  which  leads 
to  a  decrease  in  voltage  efficiency.  Fig.  7  also  indicates  that  a  higher 
methanol  concentration  may  lead  to  a  larger  current  density.  How¬ 
ever,  the  maximum  £dmfC  is  suppressed. 

By  introducing  the  V(T,  c,  j)  function  obtained  in  Section  4.1  into 
Eq.  (24),  the  Pdmfc(^  CJ)  in  the  multi-dimensional  variable  space 
was  obtained  and  shown  in  Fig.  8.  With  the  operating  temperature 


Current  density/ A  cnf  05  175  Concentration  /  M 


Fig.  6.  An  illustration  of  the  efficiency  of  a  DMFC  in  the  (T,  c,j )  space. 


5060 


Y.-J.  Chiu  et  al.  /  Journal  of  Power  Sources  196(2011)  5053-5063 


Fig.  7.  The  j-£dmfc  curves  with  various  fuel  concentrations.  The  operating  temper¬ 
ature  is  333  K,  while  the  fuel  concentrations  are:  (O)  0.75  M,  (A)  1.0  M,  (□)  1.25  M, 
and  (v)  1.5  M. 


being  333  K,  the  j-Pdmfc  curves  are  depicted  in  Fig.  9  to  illustrate 
the  influence  of  fuel  concentration  on  pdmfc.  It  is  observed  that  in 
the  region  of  c  =  0.75-1.5  M,  a  higher  fuel  concentration  leads  to 
increases  both  in  highest  j  and  maximum  pdmfc. 

Finally,  the  £dmfC(T,  CJ)  and  the  pdmfC(T,  CJ )  defined  in  Eqs.  (23) 
and  (24),  respectively,  are  illustrated  simultaneously.  Additional 
four  sets  of  j-V  and  methanol  crossover  measurement  data  were 
introduced  to  validate  the  pdmfc_£dmfc  relationship  depicted  by  the 
semi-empirical  model.  The  results  are  shown  in  Fig.  1 0  and  the  leaf¬ 
shaped  curves  are  observed.  The  symbols  in  Fig.  10  represent  the 
experimental  data,  while  the  curves  are  calculated  using  the  semi- 
empirical  model  of  Eqs.  (23)  and  (24).  A  total  of  28  points  were 
evaluated  that  shows  a  good  coincidence  between  the  predicted 
results  and  the  experimental  data.  The  relative  error  regarding  the 
Pdmfc  data  was  within  2.2%  and  that  regarding  the  £dmfc  data  was 
within  0.8%.  Accordingly,  the  proposed  models  are  verified  to  be 
capable  of  describing  the  relationship  between  Pdmfc  and  £dmfc  of  a 
DMFC  under  various  operating  conditions. 


Fig.  9.  The  j-pdmfc  curves  with  various  fuel  concentrations.  The  operating  temper¬ 
ature  is  333  K,  while  the  fuel  concentrations  are:  (O)  0.75  M,  (A)  1.0  M,  (□)  1.25  M, 
and  (v)  1.5  M. 


4.3.  Further  discussions 

One  of  the  major  objectives  of  this  paper  is  to  depict  the  per¬ 
formance  indices  (pdmfC  and  £dmfC)  and  the  operating  variables  (T, 
c,  and  j),  sin  one  diagram.  The  pdmfC  and  the  £dmfc  are  expressed 
by  explicit  functions  of  (T,  c,  and  j)  It  is  evident  that  the  operating 
range  of  current  density  varies  with  temperature  and  fuel  concen¬ 
tration.  On  the  other  hand,  the  operating  voltages  always  range 
between  zero  and  the  thermodynamically  reversible  voltage.  It  is 
more  intuitive  to  estimate  the  voltage  efficiency  using  operating 
voltage  rather  than  current  density.  Therefore,  the  operating  volt¬ 
ages  instead  of  current  densities  are  designated  on  the  pdmfc_6:dmfc 
curves  to  indicate  the  influence  of  the  electric  load  on  these  two 
performance  indices.  Whenever  a  specific  voltage  is  assigned,  the 
current  density  can  be  estimated  from  Eq.  (13)  by  adopting  the 
Newton-Raphson  algorithm  and  then  substituted  into  Eqs.  (23) 
and  (24)  to  obtain  pdmfc  and  £dmfc  for  power  density  and  efficiency 
evaluation. 


Fig.  10.  Thepdmfc  -£dmfc  relationship  calculated  from  semi-empirical  model  (curves) 
and  the  experimental  data  (symbols)  for  model  validation.  The  fuel  concentration 
is  1.5  M.  The  operating  temperatures  are:  (•)  303  K,  (a)  313  K,  (■)  323  K,  and  (▼) 
333  K. 
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Fig.  11.  The  Pdmfc -£dmfc  curves  under  various  operating  temperatures.  The  fuel  con¬ 
centration  is  1.0  M.  The  symbols  indicate  the  corresponding  operating  voltages:  (•) 
0.5  V,  (A)  0.4  V,  (□)  0.3  V,  and  (▼)  0.2  V. 

Fig.  1 1  shows  Pdmfc-^dmfc  curves  at  various  operating  temper¬ 
atures  and  a  fuel  concentration  of  1.0  M.  The  (Pdmfc  _6:dmfc)  data 
points  corresponding  to  four  levels  of  operating  voltages  are  also 
designated  on  the  curves.  Each  curve  intersects  at  the  origin  of 
Pdmfc-^dmfc  diagram.  Although  the  origin  corresponds  to  the  open- 
circuit  voltage  and  a  maximum  value  of  voltage  efficiency,  it  also 
indicates  the  current  density  is  zero;  hence  both  Pdmfc  and  £dmfc 
are  zero.  The  decrease  of  an  operating  voltage  from  open-circuit 
voltage  to  a  lower  voltage  results  in  a  decrease  in  voltage  effi¬ 
ciency,  whereas  the  current  density  is  increased  which  leads  to  an 
increase  in  the  fuel  utilization  efficiency.  Therefore,  thepdmfC-£dmfc 
diagram  shows  a  positive  tendency.  Along  a  leaf-shaped  curve  in 
the  clockwise  direction,  an  ultimate  £dmfc  is  achieved  and  followed 
by  a  decrease  in  £dmfc  because  the  decreasing  voltage  signifi¬ 
cantly  lowers  the  resulted  values  of  £dmfC  and  Pdmfc  -  Eventually,  the 
Pdmfc “^dmfc  relationship  is  back  to  the  origin  of  the  diagram  when 
operating  voltage  approaches  zero.  In  practice,  only  the  operating 
conditions  regarding  the  upper-left  portion  of  the  curve  are  consid¬ 
ered  since  they  always  achieve  a  higher  efficiency  within  a  feasible 
power  domain. 

In  Fig.  1 1 ,  the  preferable  value  of  the  operating  voltage  is  around 
0.4  V  since  bothpdmfc  and  £dmfc  approach  ultimate  values.  It  is  also 
observed  that  increasing  the  operating  temperature  is  beneficial  to 
enhance  the  power  density  of  a  DMFC  without  reducing  its  maxi¬ 
mum  £dmfc-  However,  under  a  constant  pdmfC,  increasing  operating 
temperature  could  reduce  £dmfC  due  to  the  decrease  in  fuel  uti¬ 
lization  efficiency.  Another  example  regarding  the  effect  of  fuel 
concentration  on  the  Pdmfc-^dmfc  relationship  is  shown  in  Fig.  12, 
in  which  the  operating  temperature  is  kept  at  323  K.  It  indicates  the 
maximum  £dmfc  occurs  at  about  0.4  V,  whereas  the  maximum  pdmfC 
is  within  a  range  of  0.3-0.4V.  A  higher  level  of  fuel  concentration 
causes  an  increase  in  Pdmfc-  Nonetheless,  higher  fuel  concentration 
also  leads  to  severer  fuel  crossover,  thus  significantly  suppressing 

^dmfc- 

The  remaining  run-time  tr(T,  c,  j)  is  another  index  of  interest 
for  a  fuel  cell  system.  Fig.  13  shows  the  plots  of  specific  remaining 
run-time  tr(T,  c,  j)  against  pdmfc  as  defined  in  Eq.  (26).  The  operat¬ 
ing  conditions  were  the  same  as  those  shown  in  Fig.  11.  Followed 
by  a  decrease  in  the  operating  voltage,  the  tr  shows  a  dramatic 
drop  to  a  level  corresponding  to  the  maximum  Pdmfc-  Then  it  tends 
to  approach  a  stabilized  value  that  corresponds  to  the  situation 
where  the  limiting  current  density  occurs.  It  is  noted  that  both  the 
Pdmfc “^dmfc  relationship  and  the  remaining  run-time  are  coupled 


Fig.  12.  The  pdmfc  -£dmfc  curves  with  various  fuel  concentrations.  The  operating  tem¬ 
perature  is  323  K.  The  symbols  indicate  the  corresponding  operating  voltages:  (•) 
0.5  V,  (A)  0.4V,  (□)  0.3  V,  and  (▼)  0.2  V. 


Fig.  13.  The  specific  remaining  run-time  under  various  operating  temperatures  and 
a  fuel  concentration  of  1.0  M.  The  corresponding  operating  voltages  are:  (•)  0.5  V, 
(A)  0.4V,  (□)  0.3  V,  and  (▼)  0.2  V. 

with  operating  T,  c  and  j  (or  V)  simultaneously.  When  considering 
£dmfc  and  tr,  a  fuel  cell  does  not  always  operate  at  the  condition 
that  a  maximum  power  is  delivered.  In  a  more  realistic  situation,  a 
feasible  power  demand  is  given  for  a  controlling  strategy  to  pursue 
its  optimum  value  of  efficiency.  Therefore,  how  to  determine  the 
combination  of  operating  (T,  c,  and  j)  (or  V)  is  crucial  to  a  DMFC 
system.  The  model  described  in  Eqs.  (23)  and  (24)  will  serve  as  an 
important  basis  for  this  purpose  [26]. 

5.  Conclusions 

A  semi-empirical  model  for  efficiency  evaluation  of  a  DMFC  has 
been  proposed.  The  power  density  and  the  efficiency  of  a  DMFC 
have  been  represented  by  explicit  functions  of  operating  temper¬ 
ature,  fuel  concentration,  and  current  density.  It  is  illustrated  to 
be  capable  of  depicting  performance  indices  and  operating  vari¬ 
ables  in  a  diagram.  It  provides  a  good  prediction  and  a  clear  insight 
into  these  relationships.  The  experimental  validation  of  the  pro¬ 
posed  model  is  also  verified.  Moreover,  the  remaining  run-time  of 
the  system  under  various  operating  conditions  can  be  estimated 
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Table  A.1 

The  corresponding  parameters  and  functions  concerning  Levenberg-Marquardt  algorithm. 


Unknown 

Q 

Estimated 

Corresponding 

Measured  data 

P 

parameters 

functions 

equations 

0  =  [  ri  r2]T 

2 

f=Re 

Eq.  (12) 

Area-specific  resistance 

16 

0  =  [ak  bk  Chk 

C2,k  0cr,k  ] 

5 

f=V 

Eq.  (15) 

Operating  voltage 

10 

0=[Pi  fc  ■■■ 

ft  ]T 

8 

f  =  Nmcr 

Eqs.  (16)-(22) 

Methanol  crossover  flux 

28 

Table  B.l 

A  brief  summary  of  the  fitting  parameters  defined  in  Eq.  (14). 

Parameters 

Fitting  results 

Xi(T,cf 

a-u 

Oi2i 

a3i 

olm 

Ot5i 

0t6i  ol1{ 

a8i  a9i 

GhOi 

a(T,  c) 

-4.40 

-0.25 

0.38 

-0.23 

1.19 

-0.94  0.27 

-0.48  -0.36 

0.49 

b(T,  c) 

0.15 

-0.02 

-0.31 

0.07 

-0.06 

0.59  -0.07 

0.14  -0.19 

-0.27 

Ci(T,c) 

-0.48 

0.12 

3.21 

0.09 

-1.06 

-5.61  -0.12 

0.39  0.85 

3.20 

C2(T,  c) 

88.07 

-159.70 

-137.99 

125.11 

125.60 

112.44  -35.96 

-48.79  -22.83 

-48.10 

riM,  c) 

0.59 

-0.72 

0.35 

-0.20 

2.27 

-2.10  0.38 

-1.01  -0.49 

1.52 

a  xfT,  c)  =  a-u+a2iT+a3iC+a4iT2 +a5iTc  +  a6iC2 +a7iT3 +a8iT2c  +  aQiTc2 +a^oiC3. 


by  introducing  the  proposed  model.  It  is  generally  essential  to 
practical  situations.  Since  the  model  is  stated  in  the  form  of  an 
analytical  function,  the  quantitative  prediction  of  power  density 
and  efficiency  is  straightforward  whenever  the  relevant  parame¬ 
ters  in  the  model  are  pre-determined.  Consequently,  it  can  serve  as 
an  important  basis  for  developing  controlling  strategies  of  a  DMFC 
system. 
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Appendix  A.  The  nonlinear  least-squares  algorithm 

The  well-known  Levenberg-Marquardt  algorithm  is  adopted 
for  the  fitting  of  the  unknown  parameters  in  this  paper. 
Let  &  =  [9i  02  ■■■  0q]T,  y  =  [yi  yi  yP]T,  and 

T 

f=[/i  f2  •••  fp]  denote  the  unknowns  to  be  fitted,  the 
measured  data,  and  the  estimated  results  by  the  model,  in  order. 
Based  on  an  initial  guess  of  and  p  sets  of  measurements, 
the  key  of  the  algorithm  is  to  determine  the  next  search  of  the 
unknowns  as  the  followings: 

&(»  =  ©('-d  +  [ja-ny-D  +  a(,-1)i<?]"1j('-1)T(y  -  f(,_1))  (A.i) 

where  a^\  \q  and  l  are  relaxation  factor,  identity  matrix  and  iter¬ 
ation  index,  in  order.  Furthermore,  the  Jacobian  matrix  J  is  defined 
as: 


9/i 

9/i 

9/i  1 

90, 

902 

dOq 

9/2 

9/2 

9/2 

90, 

902 

90, 

9/p 

9/p 

9/p 

_90i 

902 

90,. 

It  is  noted  that  the  finite  difference  approximation  is  introduced  for 
Eq.  (A.2)  when  the  analytical  form  of  J  cannot  be  obtained.  By  Eq. 
(A.I ),  the  unknowns  are  updated  in  each  iteration  of  the  algorithm 
until  approaches  adequately.  One  can  refer  to  [32]  for 


more  details  of  the  algorithm.  The  proposed  parameters  and  func¬ 
tions  concerning  Levenberg-Marquardt  algorithm  are  summarized 
in  Table  A.I. 


Appendix  B.  The  ordinary  least-squares  estimate 


The  unknown  parameters  o'! |,  ^2/,  • .  -,  o'10l  defined  in  Eq.(14)can 
be  estimated  by  ordinary  least-squares  method.  When  a  number 
of  (Tto  C/J  are  assigned  and  the  corresponding  parameters  xi  k  (i.e. 
ak,  Cy0  c2k,  and  rjcr<k)  are  determined  by  Levenberg-Marquardt 
algorithm,  Eq.  (14)  can  be  expressed  by  a  matrix  form: 


Xj  =  M  •  A, 

whei 
and 

M  = 


x,= 

Ki 

*i,2 

Xi.n 

\\ 

A,-  = 

[«li 

a2  i 

1  T] 

Cl 

T\ 

r,c! 

4 

n 

Tfc, 

T,c\ 

c31 

Li 

i  r2 

C2 

n 

T2C2 

4 

n 

T22C2 

Tic\ 

c3 

C2 

1  T„ 

Cn 

Tn 

TnCn 

cl 

n 

Tlcn 

TnCl 

C3 
Ln  J 

(B.l) 

J 


<*10!  \ 


(B.2) 


It  is  noted  that  n  =  1 6  in  the  illustrated  examples  of  this  article.  Then 
the  ordinary  least-squares  estimate  of  A*  can  be  expressed  as  [32]: 

Ai  =  (MTM)“1MTXi  (B.3) 


The  fitting  results  of  Eq.  (14)  are  summarized  in  Table  B.l. 
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